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ABSTRACT 

Measurements of the Faraday rotation of the cosmic microwave background 
radiation (CMBR) polarization could provide evidence for the existence of pri- 
mordial magnetic fields. The Faraday rotation could also allow the study of some 
properties of these fields. In this paper, we calculate the angular dependence of 
the Faraday rotation correlator for different assumptions about the spectral in- 
dex and correlation length of the magnetic field. We show that the helical part 
of the magnetic field does not make any contribution to the correlator. We stress 
the importance of the angular resolution of the detector in the Faraday rotation 
measure, showing that it could severely reduce the effect, even for a relatively 
large magnetic field correlation length. 

Subject headings: magnetic fields-cosmology: theory-cosmic microwave background 

1. Introduction 

Astronomical observations have revealed the presence of large-scale magnetic fields in 
the Universe. They exist in all types of galaxies (spiral, elliptical, barred, and irregular), 
in galaxy clusters, and probably in superclusters, with correlation lengths £ ~ Mpc and 
intensities B ~ yuG. The presence of magnetic fields in all gravitationally bound large- 
scale structures could suggest that they have been created after structure formation through 
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astrophysical mechanisms (as, for example, the "Biermann battery"). On the other hand, the 
detection of magnetic fields in galaxies at high redshifts, could represent a strong hint that 
magnetic fields have been generated in the early Universe (i.e before structure formation) 
by microphysics processes (for a full discussion see Grasso & Rubinstein 2001; Widrow 2003; 
Dolgov 2001; Giovannini 2004). 

If large-scale magnetic fields have a primordial origin, they could have observable ef- 
fects on the cosmic microwave background radiation (CMBR). In particular, as discussed in 
(Kosowsky & Loeb 1996), existence of a magnetic field at the last scattering surface, corre- 
sponding to a present-day value of Bq ~ 10~ 9 G, may induce measurable Faraday rotation of 
the CMBR polarization. Larger magnetic fields, at the level Bq ~ 10 _8 G, may even depo- 
larize the CMBR (Harari et al. 1997), owing to differential Faraday rotation across the last 
scattering surface. 

Clearly the measure of a non-zero Faraday rotation in the CMBR polarization would be 
extremely important, since it would support the primordial origin of cosmic magnetic fields, 
and would provide direct information about the physics of the early Universe. 

In this paper, we analyze the Faraday rotation of the CMBR polarization induced by 
a small (Bq ~ 10~ 9 G) primordial stochastic magnetic field. In particular, we calculate the 
angular dependence of the Faraday rotation measure (RM) maps for different assumptions 
about the magnetic field spectral index and correlation length. We consider the possibility 
of helical primordial magnetic fields, showing that, contrary to the suggestion of (Pogosian 
et al. 2002; Pogosian et al. 2003), the helical part of the magnetic field does not contribute 
to the Faraday RM maps at any angle. Finally, we discuss the importance of the detector 
angular resolution in the Faraday rotation measure, showing that it could severely constrain 
the actual possibility of observing the effect. 

The plan of the paper is as follows. In Section 2, we briefly review the properties 
of a stochastic homogeneous and isotropic magnetic field and we define its characteristic 
properties. In Section 3, we discuss the Faraday rotation effects and show that Faraday RM 
maps do not depend on the helical part of the magnetic field. In Section 4, we discuss the 
angular dependence of Faraday RM maps for different assumptions on the magnetic field 
spectrum. Finally, we summarize main results in Section 5. 
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2. The magnetic field spectrum 

We consider primordial stochastic magnetic field B(r/, x), created before the matter- 
radiation decoupling. Its Fourier transform, B(r/, k), is defined according to 1 

B^, k) = J rf 3 xe ,kx B I ( f | ) x), flifax) = J d 3 ke- ikx 5,(77, k). (1) 

Here 77 is the conformal time, x are comoving coordinates, and k are the comoving wavenum- 
bers. 

Under assumption of flux conservation, the magnetic field scales as a~ 2 , where a is 
the cosmological scale factor (see, e.g., Grasso & Rubinstein 2001). The magnetic field at 
arbitrary time can thus be related to its present value by B(?7, x) = B (x)/a(?7) 2 , where the 
subscript indicates today's values and we normalize the scale factor according to 0(770) = 1- 

We are interested in the effects of a statistically homogeneous and isotropic mag- 
netic field. This means that the correlation tensor of the magnetic field, C i j(r 1 ,r 2 ) = 
(B i0 (r 1 )Bj (r 2 )), is a function of r = |r x — r 2 | only and, moreover, it transforms as an 
5*0(3) tensor. In terms of the Fourier amplitudes of the field, these conditions (together 
with the fact that the magnetic field is a divergence- free field) translate into (see e.g. Monin 
& Yaglom 1975; Caprini et al. 2004) 

(B i0 (k)B j0 (V)) = ^5(k + k')[P tJ S(k) + ieijMik)], (2) 

where = 5ij — hkj, £^7 is the totally antisymmetric tensor and ki = ki/k. Here S(k) 
denotes the symmetric part and A(k) denotes the antisymmetric part of the correlator. Usu- 
ally S(k) is referred to as the magnetic power spectrum, being related to the total magnetic 
energy E B = | J d 3 xB„(x) (the volume is normalized to V — 1) by 

/•OO 

£ B (k) = 2nk 2 S(k), E B = dkE B (k). (3) 

Jo 

On the other hand, A(k) is referred to as the helical power spectrum, being related to the 
magnetic helicity H B = f d 3 x A (x) ■ B (x) by (see e.g. Vachaspati 2001): 

POO 

Ti B (k) = ATvkA(k) , H B = dkU B {k). (4) 

Jo 



^^Here and in the following we use the natural system of units in which h = c = 1. We use bold characters 
to indicate vectors, i.e. a = ~ct, while normal characters are used to indicate vector moduli, i.e. a = |a|. 
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We will assume that both S(k) and A(k) can be represented by the following simple functions: 

S(k) = S k ns e~ WK) \ A(k) = A k nA e~ {k/K) \ (5) 

which for k <C K possess a power law behavior. For large k, the spectrum is instead 
suppressed exponentially in order to have finite energy and helicity and a finite correlation 
length £, given by: 

JZ°dk(2ir/k)£ B (k) _ 2vr r[(n 5 + 2)/2] 
* - J °°dkS B (k) K T[(ns + 3)/2] ' U 

The two functions S(k) and A(k) are not completely independent, since any field configura- 
tion has to satisfy: 

S(k) > \A(k)\. (7) 

Moreover, requiring that the function (2) is analytic, one can show that the spectral index 
ris has to be even > 2, while ua has to be odd and > 2 (Durrer & Caprini 2003). 



3. Faraday Rotation 

As discussed in (Kosowsky & Loeb 1996), a primordial magnetic field could leave sig- 
nificant imprints upon the CMBR polarization through the effect of the Faraday rotation. 
Any magnetic field between the last scattering surface and the observer would rotate the 
polarization vector by the angle 2 : 

e 3 

cM> = A 2 - n e B ■ n a drj , (8) 

8n 2 mj 

where n indicates the photon propagation direction, a is the scale factor, A is the photon 
wavelength and n e is the number density of free electrons. In the assumption of flux con- 
servation, the quantity A 2 B is time independent and thus one can substitute A — > Ao and 
B — > B into equation (8). 

The Faraday rotation is proportional to the number density of free electrons, which 
evolves because of the Universe expansion and because of recombination and reionization 
phenomena. It is useful to describe these variations in terms of their contribution to the 
photon optical depth. Following (Kosowsky & Loeb 1996), we introduce the differential 



2 Here and in the following, we use the Heaviside-Lorentz electromagnetic units in which the fine structure 
constant is a = e 2 /4ir. In the previous version of this paper and in the version accepted for publication 
in ApJ, pre-factors in some formulae are incorrect. However, the final results are correct. We thank T. 
Kahniashvili for making us aware of this point. 
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optical depth f = dr/drj = n e aTd where ot = 8na 2 / 3ml is the Thomson cross section. In 
terms of this quantity, one has 

d$ = -?-T(r))\lB -ndr}. (9) 
Aire 

This expression has to be integrated along the photon path, starting from the photon last 
scatter r\ until the present time 770 ■ One obtains 

A*fa) = Ajj , (10) 

where g(rj) is given by 

q rvo 

oiv) = ^-J r dr l B o W - Vo)) n. (11) 

In the previous expression we have assumed that the observer is at the origin of our reference 
frame and n(rj' — 770) describes the photon space trajectory as a function of the conformal 
time 77. 

As a final step, one takes into account the fact that photons from a given direction 
last scattered at different times 77. One introduces, then, the visibility function 51(77) = 
f{rj) exp(— t(t))), which gives the probability that a photon observed at 770 last scattered 
within drj of a given 77, and calculates 

/to 3 rvo rvo 

d V g( V )g( V ) = —J^ d V g( V )J drj' f (77') B (n(7/ - 770)) • n . (12) 

It is straightforward to recast Eq. (12) into the form 

3 n° 

RM = — / d77W)B (n(77 , - 77o))-n, (13) 
4vre Jo 

which gives the wavelength-independent Faraday rotation measure RM in terms only of g{rj) 
and of the properties of the magnetic field B . 

In principle, by measuring the CMBR polarization at different wavelengths with a suit- 
able accuracy, one should be able to determine the Faraday RM as a function of the obser- 
vation direction and, thus, to determine the correlator 

RE! (9) = (RM(n)RM(m)) , (14) 

where n and m are two directions on the sky and cos 6 = n • rh. This correlator depends on 
the properties of the magnetic field and on the ionization history of the Universe, according 
to 

EE'(6) = (J^ J drjgir]) J drj 'g(rf) ((B (At7 n) • n)(B (A7/ m) • m)> , (15) 
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where Ai] — rj — rjo and Arj' = rf — i] . The last terms in the previous expression can be 
expressed in Fourier space as 



We remind the reader that S(k) and A(k) are, respectively, symmetric and helical terms of 
the magnetic field correlation tensor. 

In principle, in order to calculate rigorously the Faraday RM correlator RR'(8), one 
should solve the radiative transfer equations for the microwave background polarization in 
the presence of a magnetic field. The exact value of the rotation measure is, in fact, sensitive 
to the growth history of the CMBR polarization through the surface of the last scattering. 
This was done, e.g., in (Kosowsky & Loeb 1996) where the quantity RR'(0) was calculated 
in the limit of small magnetic fields (i.e. B ~ 10 _9 G), showing that, in this limit, the 
simple approach described above gives reasonably accurate results. In this paper, we use 
Eqs. (15,16), which has the advantage of simplicity and allows us to calculate analytically 
the angular dependence of RR'. This allows, in turn, to discuss in a transparent way the 
dependence of the Faraday RM maps on the magnetic field parameters, the effects of finite 
detector angular resolution and a possible role of helicity in RR'. 

In this respect, it was suggested by Pogosian et al. (2002, 2003) that the dependence 
of RR' on A(k) can be used to extract information on the helical part of the magnetic field. 
However, one can show that the helical part of the magnetic field does not contribute to RR' 
at any angle 6 3 . This is due to the fact that its contribution to the magnetic field correlator 
((B • n)(B • m)), given by 



((B • n)(B • m))hci = / d 3 k (n x m • k) A(k) exp[-ik • (nA V - riiA?/)] , (17) 



is zero for any chosen directions n and m. Indeed, by decomposing k in terms of kj_ and ky, 
where k^ and ky are perpendicular and parallel, respectively, to the plane containing the 
vectors n and rh, we get 



((B • n)(B • m))hci = / d 3 k n x m • — ) A(\k ± + ky |) exp [-ik,, • (nA^ - mArf)] . 



<(Bo-n)(B -m)) 



+ 







(18) 



3 For 8 = this follows simply from the fact that helicity affects only off-diagonal terms of the magnetic 
field autocorrelation function, as noted in (Enfilin & Vogt 2003). 
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One immediately sees that the integrand is odd in k^, and thus the integral has to vanish. 
All this shows that CMBR Faraday RM maps cannot give any information on the helicity 
of primordial magnetic fields 4 . 

In conclusion, Faraday RM maps of the CMBR only depends on the symmetric part of 
the magnetic field and we have 



RR'(0) = (J^ J d V g{ri) J drf <?(?/) ((Bo^n) ■ n)(B (Ai/m) ■ m)> 
The magnetic field correlator can be expressed as (Kolatt 1997) 
((B -n)(B -m)) symm 



symm 



n ■ m)C ± (r) + (n • £)(m • ^)(C\\(r) - C ± (r)) 



where r = hArj — mAr/, and 



2 f°° 

c ^ = wl dk£B{k) 

2 ro ° 



jo(kr) - -j 2 (kr) 



i>oo 

/ dk£ B (k)[j (kr)+j 2 (kr)], 
Jo 



3(2tt) 3 

where ji(x) are the spherical Bessel functions of the i th order. 



(19) 

(20) 

(21) 
(22) 



4. Results 

The bulk of the Faraday rotation is generated close to ?7dcc where the photon visibility 
function g(rf) is maximal. If the correlation length £ of the magnetic field is much larger 
than the thickness of last scattering surface <5?7dec ~ 10 Mpc (Spergel et al. 2003) {i.e. the 
distance travelled by a photon during the period of time in which g(i]) is sizeably different 
from zero), one can approximate the visibility function with delta function g(i]) = S(rj — r)d ec ) 
and extract ((B ■ n)(B ■ m)) symm from the two integrals in Eq. (19). As a final result, one 
obtains 

RR\e) = (^-) 2 {n ■ mC ± (r dcc ) + (n • ^)(m • ^) [C||(r dec ) - C ± (r dcc )]\ , (23) 

\ 4VTe J y ''"dec '"dec J 

where n • m = cos 9 and rd cc = (n — rh)(^dcc — Vo)- 



4 During the writing of this paper, we learned from (Kahniashvili 2004) that Kosowsky et al. are also 
studying this problem with similar conclusions. 
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For the case 9 = 0, expression (23) can be easily evaluated. We have 

M '< 0) = (2^F / * £B W = (4^ 5 » ■ (24) 



where -Bo is the average magnetic field, defined by 



B « = /(SF(2^< B »( k )' B »( k '»= (2^ B - (25 > 



This result essentially coincides with the result of (Kosowsky & Loeb 1996) and corresponds 
to the average rotation of the CMBR polarization approximately equal to 

For 9 7^ 0, the situation is slightly more complicated. After some algebra we obtain 

RR'(9) = (J^j [C ± (r dcc )cos 2 (9/2) - C||(r dec ) sin 2 (£/2) ] , (27) 

where r dec = 2(??o~ Vdcc) sin(#/2). However, in a realistic situation, the correlation length £ of 
the magnetic field is much smaller than the distance to the last scattering surface, t]q — rjdec- 
This permits us to consider the limit of small observation angles (i.e. 9 <C 1), in which the 
second term in the right-hand side of Eq. (27) becomes negligible and one obtains 

Rim = c ± (9( Vo - Vdcc )) 

RR'(0) Cj_(0) ' 1 ' 

We remark that this expression is quite natural if we consider the physical meaning of the 
function C±(r). It describes the variation of the correlator (Sj(O)Sj(r)) = C±(r) of the 
magnetic field components in a fixed direction (which in our case is the photon propagation 
direction), when r moves in a plane perpendicular to that direction (i.e. on the last scattering 
surface) . 

To obtain an explicit expression for RR'(9) one has now to specify £b(^) — 1 27ik 2 S(k) in 
Eq. (21). If S(k) can be described by the simple functions given by Eq. (5), the integration 
can be performed analytically. For the case ns = 2 we get 

where x — K(Vo ~ ? 7dec)/2. For different ns values one obtains quite similar expressions, 
of the kind exp(— x 2 9 2 )P(x 2 9 2 ), where P(x 2 9 2 ) are higher order polynomials in the x 2 9 2 
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variable [in the Appendix we present explicit expressions of C±(r) and C\\(r) for arbitrary 
values of ns]- Expression (29) can be explicitly written as a function of the magnetic field 
correlation length £, which for n s = 2 is £ = 8y/n/3K. For different ns values, one has to 
take into account that the relation between £ and K depends on the chosen spectral index, 
as it is described in Eq. (6). 

RR'(6>) 



RR'(O) 




9 (deg) 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 



Fig. 1. — Faraday rotation measure correlation RR'(9) as a function of the separation angle 
6. The three lines correspond to the magnetic field spectral index ns = 2 (solid line), ns = 4 
(dashed line) and ns = 6 (dotted line). The correlation length of the magnetic field is £ = 20 
Mpc. 

In Fig. 1 we show RR' as a function of the angle (normalized to its value at 9 = 0) for 
various choices of the spectral index ns and for a fixed correlation length £ = 20 Mpc. The 
three lines correspond to ns = 2,4, and 6. One can see that the function RR'(8) has quite 
peculiar behavior. In particular, in all cases, there are angles that correspond to negative 
values of RR'. This reflects the behavior of the magnetic field at the last scattering surface. 
Namely, it is due to the fact that the correlator {B y (0)B^(r)} ~ C±(r) (where 7 indicates the 
photons propagation direction) on the surface of the last scatter becomes negative. In other 
words, the magnetic field components £> 7 in different regions of the last scattering surface 
can be anticorrelated. 
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From the point of view of observations, it is important to note that the angular behavior 
of RR'{6) is only marginally dependent on the ng value. In particular, the point at which 
RR'(9) vanishes practically coincides in all three cases. On one hand, this indicates that the 
observation of RR'(6) could provide information on the correlation length £ that is essentially 
independent of the parameter n,s- On the other hand, it indicates that a large sensitivity is 
required to discriminate among different n,s values. 

The results presented above are obtained neglecting the last scattering surface thickness 
and assuming that the detector has a perfect angular response. These approximations are 
correct if the magnetic field correlation length £ is much larger than Srjdec and if, at the same 
time, 9% ^> a, where a is the detector angular resolution and Q$ = £/(?7o — ^dec) represents 
the angular dimension of a correlation domain on the surface of the last scatter. Both effects 
of 5rj dec 7^ and of a ^ reduce the calculated value for RR'. In this sense, Eq. (24) should 
be regarded as an upper limit for the true RR'(0) value. 

In order to estimate the role of 5?7dec in the calculation of RR'(0), one can approximate 
the behavior of g (n) around the decoupling with g(rf) = (l/v^vr^dec) exp [—(77 — r?dec) 2 /25?7^ ec ] , 
which is a Gaussian peaked at r] dcc with a characteristic width equal to Si] dcc . Introducing 
this function into Eq. (19) one can calculate explicitly the ratio between the rotation measure 
obtained taking into account the thickness of the last scattering surface, {RR'(0))$ Vdcc , and 
the "ideal" value, RR'(0), given by Eq. (24). One obtains 

(h»(o)>^._ 1 r^ exp (_^M3W> (30) 



RR(0) 2VS*Jfc ■/-«, *\ «>)L/ C,(0) 

One should note that the thickness of the last scattering surface introduces an average (along 
a line) of the function C\\(r). This is easily explained if one considers the physical meaning of 
this function. It describes (see Eq. (20)) the variations of the correlators (Sj(O)Sj(r)) of the 
magnetic field components in a fixed direction (which in our case is the photon propagation 
direction), when r moves parallel to this direction (i.e. across the last scattering surface). 
For the specific case ns = 2, we have 

CM) W 2 , s 

q$r e ~ (31) 

which gives 

/BB'frm.. r 64vr {5r ]dcc \ 2 Y 1/2 



RR(0) 



9 



For different values of the ns parameters, we obtain more complicated expressions (which 
can be calculated analytically; see the Appendix). For all the values of the ns parameter, 
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however, the ratio (RR'(0))s Vdcc /RR'(0), for 5r] dcc 3> £, always goes to zero as (5rfd.ee/0~ 1 - 
In conclusion, since <£> is proportional to (RR') 1 ^ 2 (see Eq. (26), we find that the average 
rotation angle of the CMBR photons is reduced with respect to Eq. (24) as iV" 1 / 2 , where 
N = <5dcc/£ is the number of correlation domains in the length <5?7dec- 



<RR'(0) > £ 



. 8 



. 6 



. 4 



0.2 




. 5 



1 . 5 



<r/0 € 



Fig. 2. — Effect of the detector angular resolution a on the CMBR Faraday rotation measure. 
The angle 0% represents the angular dimension of a magnetic field correlation domain on the 
surface of the last scattering. The three lines correspond to the magnetic field spectral index 
ns = 2 (solid line), ns = 4 (dashed line) and ns = 6 (dotted line). 

A similar exercise can be done to include the effect of finite detector angular resolution. 
Modelling the detector angular response with a Gaussian of angular width a (Kolb & Turner 
1990), one obtains 

(RR'(0)) a 1 / ^Kx^o-^)) 



dO 9 exp 



(33) 



RR(0) a 2 J V 2 ^V C±(0) 

where (RR'(0)) a is the rotation correlator obtained with an account of the detector angular 
resolution. We have assumed here that tr C 1. In the case ns = 2, one obtains 

-2 



(RR'(0))a 

RR(0) 



1 + 



32tt 
~9~ 



(34) 
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where the angle 0% = £/(?7o — r)a ec ) represents the angular dimension of a correlation domain 
on the surface of the last scattering. For other ns values, we obtain different expressions (see 
the Appendix) which, for a ^> 0^, always go to zero as (9^/a)^ 4 . Thus, the rotation angle 
$ is reduced because of effect of the angular resolution as (a/9^)~ 2 . We remark that this 
behavior could introduce severe limitations to the actual possibility of observing the Faraday 
rotation of the CMBR polarization. In particular, in Fig.2, we show [(RR'(Q)) a / RR^Q)] 1 / 2 
as a function of the ratio a/0% for different n s values. One sees that even for small values of 
a/0£ the reduction of the effect can be substantial. 

We point out, finally, that in contrast to the previous case, $ is reduced as iV" 1 (and not 
as N a 1 ^ 2 ), where N a = (a/9^) 2 is the number of correlation domains in a (two-dimensional) 
region at the last scattering surface of angular dimension a. This is due to the peculiar 
behavior of the Cj_ functions or, in other words, to the fact that magnetic field components 
in different regions of the last scattering surface can be anticorrelated. 

5. Conclusions 

In this paper we have analyzed the Faraday rotation of the CMBR polarization induced 
by a primordial stochastic magnetic field. In particular, we have calculated the correlation 
between the Faraday rotation measures, RR'(9), as a function of the separation angle between 
observation directions, 0, for different assumptions about the magnetic field spectral index 
and correlation length. Here we summarize here our main results. 

i) The helical part of the magnetic field does not contribute to RR'{9) at any angle 9. 
This means that Faraday RM maps can provide information only on the symmetric part of 
the magnetic field. 

ii) In the approach described in Sec. 3, neglecting the last scattering surface thickness 
and the detector angular resolution, the Faraday RM maps RR'(9) can be calculated an- 
alytically. We have provided an analytic expression for RR'{9) for arbitrary values of the 
magnetic field spectral index and correlation length. 

Hi) The last scattering surface scale thickness 5rj dec reduces the Faraday rotation angle 
with respect to the simple estimate given by Eq. (26). In the limit 5r}d ec £, where £ 
is the magnetic field correlation length, the rotation angle is reduced as (5^dcc/£) -1 ^ 2 (see 
Eq. (32)). 

iv) The detector angular resolution a could drastically reduce the possibility of observing 
the Faraday rotation. In the limit a 9%, where 0% is the angular dimension of a magnetic 



-13- 



field correlation domain, the rotation angle is reduced with respect to Eq. (26) as (cr/9^) 2 
(see Eq. (34)). As shown in Fig. 2, the reduction can be substantial even for small values of 



We would like to thank D. Comelli for helpful discussions. 
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A. The C±(r) and C\\(r) correlators 

In this appendix, we calculate the correlators C±(r) and C\\(r) defined by Eqs. (21) and 
(22), for the functional form of the magnetic power spectrum given by Eq. (5). To this end, 
we introduce the function 

_ f °° dfcfc".+V«* a /* a [jojkr) - j 2 (fcr)/2] 

related to Cj_(r) by 

W0) = |||. (A2) 
It is straightforward to obtain the recursion formula for L ns (a) 

Qra 

L ns {a) = A m —L 2 (a) , (A3) 



where the "generating function" L 2 (a) is given by 



-K 2 r 2 /a / K, 2 r 2 

a 0/z \ a 



L ^) = —^l l - — )i (A4) 



and 

n s , # 3(-l) n 2 n 

m = — - 1, re = — , A n = ■ (A5) 

2 2' (2n + 3)!! v ; 

(we remember that n$ > 2 is an even natural number). 

Taking into account Eqs. (A3)-(A4) we cast Eq. (A2) in the form 

^44 = e^P^reV) , (A6) 
Cj_(0) 2 v ; 

where P n (x) turn out to be polynomials of (n + l)st degree defined by 

P n (x)=A n e x ^^-(l--)\ a=1 , n = 0,l,2,... (A7) 



c7n" n"' 2 V (\ 
The first three P n polynomials are 



P (x) = 1 - x, 

Pi (:r) = 1-^x + lx 2 , (A8) 
5 5 

. . 13 8 „ 4 , 
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Following the same procedure, we find for C\\ (r) the expression 

C{r) e-^Q^V), (A9) 



C||(0) 

where Q n (x) are polynomials of nth degree given by 

Qn e ~x/a 

Q n (x)=A n e x —-^ JY \ a=1 , n = 0,l,2,..., (A10) 

where k and A n are the same as in Eq. (A5). We give the expressions for the first three Q n 
polynomials: 

Qo(x) = 1, 

Qi(x) = l-\x, (All) 



4 4 
Q 2 (x) = 1 - - x+ — x 2 . 

Now, having the expressions for C±(r) and C\\(r), we can easily calculate the correlators 
defined by Eqs. (30) and (34). Taking into account Eqs. (A9)-(A10), we cast Eq. (30) in the 
form 

(RR'(0))s Vicc _ A d m 1 

RR'(0) m da™ a V« + (^dec) 2 

while, inserting Eqs. (A6)-(A7) into Eq. (34), we get 

(i^(0)) d- 1 

RR'(0) ~ Am da- + 2 X V2) 2 |q=1 ' (Ai ' j 

where we remember that x — K(Vo — ?7dcc)/2, and m is given by Eq. (A5). 
Finally, we give the asymptotic expressions 

(RR'(0))s Vdcc __ 3 /<W~ ] 



^^IfJ ' (A14) 

where we remember that £ is the correlation length (see Eq. (6)), the coefficients C m are 

= 3 r(m + l/2)[r (m + 5/2)l a 

(2tt) 4 [r(m + 2)] 4 ^ ' 

and ^ = £/(r/ — r/ dcc ) is the angular dimension of a correlation domain on the surface of the 
last scatter. 
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